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neutral gauge bosons are reconsidered when the scale of the 3-3-1-1 breaking is much larger 
than that of the ordinary SU(3)c <8> SU(3)l <8> U(l)x (3-3-1) breaking. We investigate how 
the 3-3-1-1 model generates an inflation by identifying the scalar field that spontaneously 
breaks the U( 1)n symmetry to inflaton as well as including radiative corrections for the 
inflaton potential. We figure out the parameter spaces appeared in the inflaton potential 
that satisfy the conditions for an inflation model and obtain the inflaton mass an order of 
10 13 GeV. The inflaton can dominantly decay into a pair of light Higgs bosons or a pair of 
heavy Majorana neutrinos which lead, respectively, to a reheating temperature of 10 9 GeV 
order appropriate to a thermal leptogenesis scenario or to a reduced reheating temperature 
corresponding to a non-thermal leptogenesis scenario. We calculate the lepton asymmetry 
which yields baryon asymmetry successfully for both the thermal and non-thermal cases. 

PACS numbers: 12.60.-i, 98.80.Cq, 98.80.Ft 


‘Electronic address: 
' Electronic address: 
* Electronic address: 
§ Electronic address: 


dthuong@iop .vast. ac. vn 
pvdong@iop.vast.ac.vn 
cskim@yonsei.ac.kr 
ntthuy@iop.vast.ac.vn 



2 


I. INTRODUCTION 

Cosmological inflation is a popular postulate for the early universe. It can solve the difficulties 
of the hot Big Bang theory and provide the predictions for quantum fluctuations in the inflating 
background. In order to recover the conventions of the hot Big Bang theory and to know how 
the universe is reheated, we must understand what is the inflaton field, and how it is connected to 
particle physics. These problems were first investigated with the chaotic inflation scenario by Linde 
[TJ . According to this scenario, the inflation may begin even there was no thermal equilibrium in 
the early universe. It can occur in a theory with a very simple potential such as V{4>) oc cj) 2 . There 
is no limit to the theory with a polynomial potential: Chaotic inflation occurs in any theory where 
the potential has a sufficiently flat region [T|. On the other hand, the recent measurements of 
B modes by BICEP2 collaboration [2] have yielded very interesting results, which could be the 
direct measurements of quantum gravitation excitations from the early universe. The ratio of the 
tensor and scalar is measured as O.I 6 +O 05 . Combined with the Planck and WMAP measurements 
suggests that the inflation model must be a larger field model. Hence, the inflationary scenario 
does not work on the framework of the Standard Model (SM) without a non-minimal coupling to 
gravity. 

Furthermore, what is the origin of matter-antimatter asymmetry in the universe? The neutrino 
experiments such as Super-Kamiokande [3j, KamLAND [3] and SNO [5] have confirmed that the 
neutrinos have small masses and large flavor mixing. According to the Planck mission team, and 
based on the standard model of cosmology, there exits dark matter (DM) which lies beyond the 
SM. All the experiments call for extensions beyond the SM. One way to extend the SM is to 
expand the gauge symmetry group. There exists a simple extension of the SM gauge group to 
SU(3)c <8> SU(3)l <8> U(l) x , the so-called 3-3-1 models. These models can explain the following 
issues paj; 

• Why the electric charges are quantized? 

• Why there are only three observed families of fermions? 

• Why top quark is oddly heavy? 

• Why the strong CP nonconservation is disappeared? 

• 3-3-1 models can provide the neutrino small masses as well as candidates for the DM HHS]. 
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There have recently emerged an extension of the 3-3-1 models, based on the SU (3)c ® SU (3) 
U(l)x <8> U( 1)jv (3-3-1-1) gauge group, which not only contains all the good features of the 3-3-1 
models as mentioned M, but also has the following advantages; 

• The B — L number is naturally gauged by combination of the SU(3)l and U(1)n charges. 
It leads to an unification of the electroweak and B — L interactions. 

• The right-handed neutrinos appear in the model as fundamental fermions that solve the 
small masses of neutrinos through a type I seesaw mechanism. 

• There exists a W-parity symmetry as a (Z 2 ) remnant subgroup of the gauge symmetry. 
Almost all the new particles have wrong lepton numbers transforming as odd fields under 
W-parity. The lightest wrong lepton number particle is identified to the DM. Because of 
W-parity conservation, the model can work better under experimental constraints than the 
3-3-1 models. 


Other highlights of the 3-3-1-1 model is that the energy scale of the symmetry breaking U(1)n, 
which can happen at a very high scale like the GUT one [ 8 j. The inflationary scenario can be linked 
to U(1)n breaking and driven by the Higgs < j> potential. Due to the local gauge U(1)n symmetry, a 
radiative correction to the inflaton potential can arise from the coupling of inflaton with the U( 1)jv 
gauge boson (Z 2 ). There exist the couplings of inflaton 4> with right handed neutrinos and Higgs 
triplets, which also contribute to the inflation potential. We would like to stress that the well- 
known advantages of a spontaneously broken gauge U( l)v symmetry include a seesaw mechanism 
for the neutrino physics [ 8 j. The presence of the right-handed neutrinos that directly interact to 
the inflaton may be compatible with the leptogenesis scenario. The aim of this work is to show that 
the chaotic inflationary scenario can be driven by the singlet Higgs potential. We also focus on 
the leptogenesis happened after the inflation through the couplings of the right-handed neutrinos. 

Our paper is organized as follows: In section[TTJ we briefly review the 3-3-1-1 model and specially 


concentrate on the Higgs and gauge boson spectra in the large A limit. In section III we present 
the inflation model by assuming the singlet Higgs </> as an inflation field. The leptogenesis related 
to the matter-antimatter asymmetry of the universe and neutrino properties is studied in section 
IV Finally we summarize our works in section [VJ 
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II. BRIEF DESCRIPTION OF THE 3-3-1-1 MODEL 


The fermion content of the 3-3-1-1 model is given as mm 


YaL = 


e a L 

V (NaRY ) 


QaL — 


daL 
UaL 

\ D aL J 


-2/3), 


-1,-1), 

^ U3L ' 

Qzl = 

d 3 L 


\ U J 


(3,3,1/3, 2/3) 


u aR ~ (3,1, 2/3,1/3), d aR ~ (3,1, -1/3,1/3), 
U R ~ (3,1, 2/3, 4/3), ~ (3,1, -1/3, -2/3), 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 


where the quantum numbers in the parentheses are defined upon the gauge symmetries ( SU(3)c , 
SU(3)l, U( 1)x , U(1)n), respectively. The family indices are a = 1,2,3 and a = 1,2. N R , U and 
D are the exotic fermions, which have incorrect lepton numbers. The other fermions have ordinary 
lepton numbers. Note that the neutral fermions N R are truly sterile since they do not have any 
gauge interaction, which contradicts to the v R ones as usually considered. 

To break the gauge symmetry, one uses the following scalar multiplets [8|: 



( pt\ 


^ Vi ^ 

p = 

too 

~ (1,3, 2/3,1/3), ?7 = 

% 


X / 


\ ^3 / 


xp 



X = 

X2 

~ (1,3, -1/3, -2/3), 

~ (1, 


U§ J 




(1,3,-1/3,1/3), 


( 6 ) 


(7) 


with the VEVs that conserve electric charge and i?-parity being respectively given by 


(p) = ^(M:°) T > Yl) = ^{u, 0,0) t , 


(x) = ^=(0,0,u) T , {</>) = -J= A. 


The pattern of the symmetry breaking of the model is given by the following scheme 
3-3-1-1 ix M v) SU(3) C ® U(1) q ®U(1) b _ l ^SU(3) c ® U(1) Q ®P 
where the electric charge Q, B — L and matter parity P take the forms, 

- 1 t h + x, b-l = -- 2 =t 8 + n, P = (- 

v3 v 3 


( 8 ) 


(9) 


Q = t 3 


( 10 ) 
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Here, Tj (i = 1, 2,3,8), X and N are the SU(3)l, U( \)x and U(l)jy charges, respectively. 

The Lagrangian of the 3-3-1-1 model is given by [8j: 

C = Y. Vi'fDp* + Y 

fermion multiplets scalar multiplets 

--GiuuG^ - 1 A iuv A? v - X B UV B^ - 1 C UV C> W 

^ L H‘ U l ^ L h LU l ^ h LL/ ^ h LL/ 

~V{p,r},x,(/>) + ^ Yukawa) (ii) 

where the Yukawa Lagrangian and scalar potential are obtained Eli as follows 

^Yukawa = Kb^aLP^R + Kb^aLWbR + h'^aR^bR^ + h U Q 3 LXUr + h®pQ aL X* D p R 

+KQ 3 LVU a R + h d a Q 3 L pd a R + h d aa Q aL p* d a R + h^ a Q aL p*u a R + H.C, ( 12 ) 
V{p,V,X,<t>) = plp^p +pWx +pWv +Mp^p ) 2 + Mx^x) 2 + hi^r]) 2 
+A 4 (/? t p)(x t x) + A b{p ] p){rj ]/ n) + Ae(x t x)(?7 t ?7) 

+A7 (pYXxV) + A s{p ] p){p ] p) + Ag(x t r?)(?? t x) + (. f^ mnp VmPnXp + H.c.) 
+hW + A(^V) 2 + A10 4 >){p ] p) + An(<?!> t ^)(x t x) + A 12(^1^) (??!??)• ( 13 ) 

Because of the 3-3-1-1 gauge symmetry, the Yukawa Lagangian and scalar potential as given take 
the standard forms which contain no lepton-number violating interactions. 

The fermion masses that result from the Yukawa Lagrangian have been presented in [8j. The 
phenomenology of the 3-3-1-1 model with the A scale of the U( 1)n breaking comparable to the u 
scale of the 3-3-1 symmetry breaking has been studied in [9j. Below, we will compute the physical 
states and masses for the scalar and gauge sectors in the limit A u, which is needed for our 
further analysis. 


A. Scalar sector 


In this part, we identify the physical particles in the scalar sector. We expand the neutral 
scalars around their VEVs 0 such as 

/ nf \ f*&(u + S 1 +iA 1 )\ ( 4 jRl±i.A'X \ 


P = 


Pi 

(v + S 2 + 1 A 2 ) 

pi 


\v = 


% 


;x = 


V :/ 2 ( S !i + iA :i) ) 


X.2 

^ + £>3 + ^ 3 ) J 


; (14) 


= (A + S*4 + iA±). 
\/2 


rsj 


(15) 
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In the scalar sector, all scalar fields with W-parity even, Si, S 2 , S 3 , S 4 , mix via the mass matrix 
such as 


Mf = 


( 2X 3 U 2 - X 5 uv+^fu Xeuu + ^fv X 12 uA \ 

A 5 UV + 4 =fu 2 X 1 V 2 - -4/4r A A uv + 4= fu X 10 vA 


Xqiiuj + A=fv X^ujv + A=fu 2X 2 uj 2 - A n wA 


V 


Ai 2 uA 


AiquA 


AnwA 


2XA 2 


(16) 


We assume that A u ~ —/ S> u,v then the mass matrix given in Eq. (16) has form as 


where 



(17) 


A = 2XA 2 , 


(18) 


B — ^ A 12 UA AiquA AnwA ^ 


(19) 


C = 


(W-ig X,uv + % A) 
A5 UV + 2X\v 2 - X A vu + 


\/2 

fv 


V 2 v 

fu 


a/2 

fuv 


( 20 ) 


V + $ A 4 ucu + $ 2 W - Jg y 

Since (A 3> — /, w > M,r), we get A^> B,C. The matrix given in © can be diagonalized by 
using block diagonalizing method. The approximately unitary matrix U, 


( 


U = 



V 


1 

0 

0 

u\l 2 
' 2AA 


0 

1 

0 

' 2AA 


0 

0 

1 

' 2AA 


mAj2 \ 
2AA 

^Aio 

2AA 

2AA 

1 


( 21 ) 


transform Mj into approximately block-diagonal form: 


u ] m 2 s u 



( 22 ) 


In the limit (A 3> —f,uj S> u,v), U ~ I then H 3 ~ S 4 gets mass m ^ 3 = 2AA 2 . Si,S 2 ,S 3 are 
mixing with the mixing mass matrix obtained as 
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C-B^A~ X B = 


2A 3 « 2 - ft - ^ ^uv + f 2 -si^ A 6 W + A _ S^4n \ 


A 5 uu + - M,,A 2 1 ° Al2 2 Aiw 2 - - ^Aio A 4 uw + - ^AioAii 


2 A 

y XqULO + A 4 UW + 


2 A 

Aiq/) 

2 A 


f u _ w*jAiqAu 

\/2 


2 A 


2 w-^-^y 


2 A?, 


(23) 


At the leading order (—f,u S> u, v), the mass matrix given in Eq. (23 ) can be rewritten as 

\ 


/ f VUJ f u 

y/^U y /2 

fu fuu 

y/2 s/2v 

0 2A 2 cu 2 - ^ 


0 


0 


V 0 


(24) 


2 A 


The physical fields with respective masses can be written as 


uS 1 + vS 2 2 _ 

B — / 9 , 9 ~; m H ~ 

VM + V z 

-vSi+uS 2 2 f(u 2 + v 2 )u 

B 1 — 7 n = ITT'Hi ~ ~ 


\Ju 2 + V 


H 2 = S 3 , m Ho = 


2 ’ ~ ni y/ 2 uv 

(4AA 2 - A^)w 2 


2A 


(25) 


In the new basics, (H, Hi, H 2 ), the squared mass matrix given in (23) can be written as 


where 



(26) 


a = 


u 4 (4AAi — A 2 0 ) — u 4 (A 2 2 — 4 AA 3 ) — 2u 2 u^(AioAh — 2 AA 5 ) 


2 (u 2 + u 2 )A 


(27) 


mu(h 2 (Aio(—A io+An)+A( 4 Ai — 2A5 ))+i/ 2 (—A10 Aii+A^—4AA3+2AA5)) 

j^r _ | 2{u 2 -\-v 2 )X 

2y/2fuvX— cd(ii 2 (AioAii—2AA4)+ii 2 (An A 12 —2AAe)) 
2y/u 2 -\-v 2 X 


(28) 


-y/2f (ii 2 +i; 2 ) 2 a;A+n 3 t; 3 (—(Aiq—A i 2 ) 2 + 4 A(Ai+A 3 —A 5 )) y/2 f {u 2 —v‘ 2 )X-\-uvuj(— A 10 An+An A 12 + 2 AA 4 —2AAg) 


_ | 2uv(u 2 -\-v 2 )X 

\/2/(ti 2 — v 2 )X-\-uvuj(— A 10 An + An A 12 + 2 AA 4 —2AAe) 
2\Ju 2 -\-v 2 X 


2y/u 2 -\-v 2 X 

\/2fuvX+uj 3 (A 2 2 — 4 AA 2 ) 
2cjA 


(29) 


Since —/, w u,v, we get A' B',C'. If we kept explicitly the O(^), the H\,H 2 ,H Higgs 
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bosons can gain mass by using block diagonalizing method as 


m' 2 

m H 1 

m' 2 

m H 2 


m' 2 

m H 


9 „ 'U, V , 

m H ] + ), 

UJ 

m H 2 + °(~)’ 

UJ 

v 4 ( 4 AAi - Af 0 ) — tt 4 (Af 2 — 4AA3) — 2u 2 n 2 (AioAi2 - 2AA5) 

2 {u 2 + v 2 )X 

+ 2V2(u> + Wi, - 4 AA 2 ) (m ° + mi i + m2 &' 


(30) 


where 


m 0 


= V / 2 (c 2 (AioAh — 2AA4) + u 2 (AnAi2 — 2AAg)) 2 , 


mi — 8rmA(n 2 (—AioAn + 2AA4) + u^(—A11A12 + 2AAg)), 
m2 = 8v / 2w 2 n 2 A 2 . 


(31) 


For the remaining fields in the pseudoscalar sector, the mass spectrum is similar to that of work 
given in [9]. Let us give a brief result. 

• The pseudoscalar A 4 is massless and is identified to the Goldstone boson of Z^. 

• Two other fields are massless that are identified to the Goldstone bosons of Z and Z' 

—uA\ + VA 2 


G z = 


y/u 2 + • 


G ' = V 1 Ai+v 1 A 2 ) + (u 2 + v 2 )A 3 

y/ ( U~ 2 + V~ 2 + U J~ 2 )(u ~ 2 + V~ 2 ) 


One neutral complex Goldstone boson , Gx = , that is eaten by X gauge boson. 

One neutral complex Higgs , namely H' = U J 1 ^' I 2 with the squared mass m 2 H , = (|Ag — 
-f^){u 2 +u 2 ). 

y/zlLOJ 7 v 7 

One physical pseudoscalar (A) with mass 

/ u 2 v 2 + u 2 u > 2 + v 2 oj 2 


m\ = 


%/2 


UVLU 


(33) 


and the physical state respectively 


A = 


u L A\ + v l A 2 + uj l A 3 


y/u 2 + v 2 + c jJ 2 

For charged scalars, the mass spectrum is seminar to that of work given in [8j. 


(34) 


H = vx 2 + up 3 


y/v 2 + 


OJ 


2 ’ 


H7 = 


yy -2 + u Pi 

y/u 2 A v 2 ’ 


(35) 
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with respective masses 


m fu = ( ^ a 7 - 


-4^-) (v 2 +u 2 ), m 2 H =(^X 8 -^^)(u 2 + v 2 ). (36) 

V2uu;/ \2 \J 2 uv, 


The model contains two massive charged Higgs and two massless Higgs that are identified 
to the Goldstone bosons of Y and W bosons. 


GZ = 


ux 2 - v P:s 


\Jv 2 + UJ 2 


n- _ urj 2 - cpi 

^ _ r~ 2 ~\—2 
v it + v z 


(37) 


B. Gauge sector 


In this section, let us consider the gauge boson spectrum. The mass Lagrangian is given as 


Ggaugemass = ( 0 ) 0 ) ^^) {d A a /j,T a — —gxBg — —gxCg) ( 0 , 0 , —^=) 

+(^|) 0, 0)(gA afJj T a — -gxBg + - gNCg) 2 {-j =, 0,0) T 
+(0, ) 0) (gA ail T a + -gxBg + -gNCg) 2 (0, —j=, 0) T 


+2 {gNC^A) . 


(38) 


Let us denote the following combinations 


w ± ~F 

vv g 


V2 


ft = 

/t 


^6/i T 

V2 ' 


The non-Hermitian gauge bosons , Y)^ have the following masses 


Alw = \d 2 (u 2 + u 2 ), Ml- = p 2 (v 2 + uj 2 ). 


(39) 


(40) 


It is worth noting that A^ and A§g gain the same mass. Therefore, these vectors can be combined 
the following physical states 

A^f! — %A^g 


X° = 


V2 


(41) 


and its mass is given: 


M 2 x = ^g 2 (u 2 +uj 2 ). 


(42) 


There is a mixing among A$g, A$g, Bg, Cg components. In the basis of these elements, the mass 
matrix denoted by M 2 is given as follows 
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( \{u 2 + v 2 ) 


9 


2 V 3 


2 V 3 

Uu 2 + v 2 + 4<u 2 ) 


tl(u 2 +2i> 2 ) 

3 

t\(u 2 —2(v 2 +oj 2 )) 

3vl 


t 2 {u 2 —V 2 ) 

3 

t 2 (u 2 +v 2 + 4w 2 ) 
3V3 


2 _ t l{ u 2 -2^+u 2 )) 2 t 2. {u 2 +Av 2 +u 2 ) _ 2^ („2 _ 2(t; 2 + w 2 )} 

t2(u 2 — V 2 ) t2{u 2 +V 2 -\-4uJ 2 ) 

3 375 


V 


- 2(v 2 + w 2 )) ^(w 2 + V 2 + 4(w 2 + 9A)) J 


, (43) 


where t\ = g x /g, h = gN/g- 


The mass matrix in ( |43[ ) contains one exact zero eigenvalue with the corresponding eigenstate 
as follows 


A " ~ v'sTiif ( M3 '‘ v?' 48 " + B ") ' <44) 

It is worth to notice that A^ is the combination of Az^^A^, and B^ without contribution of the 
new gauge boson The factor t\ can be expressed in term of the sine of the weak mixing angle 
s\v by identifying the coefficient of the eey vertex with the electromagnetic coupling constant e, 
similarly as the analysis in m- We get 


1 1 = 


V%sw 

3 — 4 Syy 


(45) 


The diagonalization of the mass matrix is done via three steps. In the first step, in the base of 
A^, Z^. Z' jn C tI , the two remaining Z^, Z' gauge vectors are given by 

7 \/3 + 1 2 ti(y/3tiAsfj, — 31 ?^) 

Z,i — — , ,, + 


M V / 3 + 4 tj 

V3 


Z'„ = 


^ 8 u + 


In this basis, the mass matrix M 2 becomes 


V / 3T^f^/3T4^ 2 ’ 

h B„. 


(46) 



(47) 


where M ' 2 is the 3x3 mixing mass matrix of Za , Z' C^ gauge bosons given as 


/ 


(3+At 2 )(u 2 +v 2 ) 


x /3+4t 2 ((-3+2t 2 )« 2 +(3+4t 2 > 2 ) 

6(3+t 2 ) 


\J 3+4t 2 t 2 (u 2 —v 2 ) 
3^/3+t 2 


2(3+tf)^ 

■\/3+4t 2 ((—3+2t 2 )u 2 +(3+4t 2 )v 2 ) (3-2tf ) 2 u 2 +(3+4tf ) 2 u 2 +4(3+t? ) 2 ^ 2 i 2 ((3-2tf)if 2 +(3+47j> 2 +4(3+t? )u 2 ) 

6(3 +ff) 18(3+t 2 ) 9y/3+t 2 


y/3+4 t\t 2 {u 2 —v 2 ) 
3y/3+t 2 


t 2 «3-2t 1 )u- + (3+^ +M3 +h )u ) 2 t 2^ u 2 +v 2 + ^ UJ 2 + 9A 2 )) 


( 48 ) 
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The matrix given in (48) can be diagonalized by using block diagonalizing method. In new basis 


(Z^. Z'^, ZN), the mass mixing matrix is given as 

' A 0 
0 m 2 zN 

where A is the 2x2 matrix 

(3+4t 2 )(u 2 +v 2 ) (n( v 4 \ 

2(3+tl 2 ) 

V^3+4Tj((— 3+2tj)«“+(3+4fj)tJ 2 ) | fn ^v 2 U} 2 \ 

6(3 +tf) +t2(^j 


A= 9 - 


V3+4t 2 ((-3+2i 2 )« 2 +(3+4t 2 > 2 ) | /nr 
6(3+t 2 ) + 


A 2 


(3-2t 2 ) 2 iA+(3+4t 2 ) 2 u 2 +4(3+t 2 ) 2 oj 2 

18(3+i 2 ) 


+ 0 ( 


A?) 


(49) 


(50) 


m 


U ^ ±g 2 t 2 A 2 . 


(51) 


The new basis (Z^,Z^, Zjj) is related to the basis (Z^, Z'^, C^) as following 


where 




1 1 0 -ei 



% 

= 

0 1 —e 2 





(h £2 1 / 


UJ 


ei = 


3^3 + At\(u 2 — v 2 


C2 = 


2-^/3 + t\t 2 (u 2 + v 2 + 4(cu 2 + 9A 2 )) ’ 
(3 — 2t\)u 2 + (3 + 4i 2 )u 2 + 4(3 + t 2 )oj 2 
2-^/3 + t\t 2 {u 2 + v 2 + 4(cj 2 + 9A 2 )) 


In the limit A w u, v, 


Za ~ z„, z’ ~ z' z" ~ c„. 


(52) 


(53) 


(54) 


The new heavy gauge boson Z^ is imbedded to the gauge group U( 1)n- It approximately does 
not mix to other gauge bosons. Z p and Z' fl are mixing of the two physical field Z 4 , Z 2 . 


z l = cos^-sinf.2^, Z^sin^ + cos^, 
u 2 + 4u 2 + a; 2 


2 9 ( 2,2, 

m z ,^ J lu+u+ 3 _ 4s2 


w 


- 4 


yJcyyU 4 + V 4 - C2WV 2 00 2 + C^W 4 + U 2 {-C2WV 2 + (—1 + 2s^ / )cj 2 ) 


(3 - 4 S 2 ) 


.r ( ?/2 , , ,2 . W 2 + 4u 2 + W 2 
8 


m z2 - (u +u + 3 _ 4s2 




+ 4 


•y/c^u 4 + v 4 - C2iyf 2 w 2 + C^CJ 4 + U 2 (—C2WV 2 + (—1 + 2s‘\ v )uj 2 ) 


(3 - 4s 2 ) 


(55) 
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whPT-P tan - ^3-4 s 2 w (c 2W u 2 -v 2 ) 

wheie tan_ s - (( _i +2 ^ )u 2 _ Cw , J 2 +2c 4 ^ a; 2 ) ■ 

If we assume ui S> u,v, then tan 2£ —> 0. We get 


i 7 2 g 2 (u 2 + v 2 ) 

- Zp, m z 1 ~ 


Z Ji ~ 


Z , 2 ~ 


Z 7 


4 r 2 

n 2 2 , .2 
2 y c vk u; 
m Z 2 — 


(3 4 Sjy ) 


(56) 


The gauge boson Z* is identified as Z„ in the standard model. 


III. GENERATION OF INFLATION IN THE 3-3-1-1 MODEL 


We would like to note that the scalar singlet cj> is completely breaking U{ 1)n- The vacuum 
expectation value (VEV) < <fi > can stay at the same scale as w’s scale and the interesting phe¬ 
nomenology of the model at TeV scale was studied in [9]. In a different situation, this VEV can 
be very high that can be integrated out from the low energy effective potential and a new gauge 
boson Zn decoupling from the gauge boson spectrum. In this part, we expect that the VEV of <fi 
is very high and consider the singlet scalar </> plays the role of inflaton field. The potential for (j> at 


the tree level can be read off from Eq. (13) as 


y<t> = + A(0V) 2 + Aio + An (0 t 0)(x t x) + A 12 


(57) 


Due to the larger VEV of cf>, the interaction terms of the singlet scalar Higgs and the ordinary 
3-3-1 model Higgs triplets can be ignored. During inflation, we get 


Vcj, = g 2 <ft(t> + \((^(f)) 2 . 


(58) 


This potential is taken part in the chaotic inflation. However, the inflaton field has coupling to 
the matter fields which allow it to make the transition to hot bing bang cosmology at the end of 
inflation, namely 




(59) 


We take into account quantum corrections to following the analysis of Coleman and Weinberg 


mi 




A 2 


( 60 ) 
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where i = v aR , </>, C' /t . x, p, V- 


m 


VaR 


= -2 


= 2 Ai 0 $ 2 ; 


"4 = 2(4 + 3A<h 2 ); = 8^$ 2 ; 

m 2 = 2A n 4> 2 ; m 2 = 2Ai 2 4> 2 . 


We get 


Kff = 


1 <j > 2 

j~2{[ *^2 ^(4(() 4 + 192^ + 4(A 2 0 + A 2 ] + A? 2 )]$ 4 hi — ^ 


647r 2 

l 

, O r . \ 9 \ 9 ! “t - 3A<£> 2 

+4(/i 2 + 3A$ 2 ) 2 In A2 -} 

= ^{«* 4 1 + 4 (4 + 3 A* 2 ) 2 In }, 


A 2 


where 


(61) 


(62) 


a — 2[—32 ^(h'i) 4 + 192<7^ + 4(A 2 0 + A 2 X + A 2 2 )]. (63) 

i 

We identify the inflaton with the real part of the B — L Higgs field, <h = \/21Z[(j)\. In the leading- 
log approximation, we obtain 

F(4>) = Vtree + E eff ~ ^4> 2 + ^ + E eff . (64) 


We would like to remain that the inflation occurs as the inflaton slowly rolls to the minimal 
potential. The inflationary slow roll parameters are given [T2] by 


1 / t/7 \ 2 


fV' 

e ($) = - mj, ly 


'V" 

V($) = m P ( 


c 2 ($) = 4 


V’V" 

“h 2 " 


(65) 


where mp = 2.4 x 10 18 GeV and a prime is denoted as a derivative of T. The slow roll condition 
means that e(4>) <C 1, | rj($) |<C 1,£(4>) <C 1. In this limit, the spectral index n s , the tensor to 
scalar ratio r (a canonical measure of gravity wave from inflation) and the running index a can be 
written as 


n s = 1 — 6e + 2rj, r = 16e, a = IGep — 14e 2 — 2£ 2 . 


( 66 ) 


The spectrum index n s is estimated by BICEP2 experiment [2], Planck |13| and WMAP9 [T3j 
measurements. It is closed to 0.96. The tensor to scalar ratio is proven by BICEP2 m , r — 0 .20t°;° 5 7 
while the Planck and WMAP9 experiments gave the bound r < 0.11(0.12). 

The number of e-folds is given by 


N = 


f$0 


Vd<5> 

~V r ' 


(67) 
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where 4> e is the inflaton value at the end of inflation and defined by max(e(<h), r/(<I>), £($)) = 1. 4>o 
is the inflation value at the horizon exit. The value of N is around 50 — 60 and depends on the 
energy scale during inflation. 

The amplitude of the curvature perturbation is given as follows 


Ai = 


V 


n 247r 2 mf,e(<4*) ’ 


( 68 ) 


-9 


The value of curvature perturbation should satisfy the Planck measurement [15]: = 2.215 x 10 

at the scale k o = 0.05 Mpc~ l . 

Let us study parameter space of /i,A,A,a appeared in the potential V(3>). If /r 2 3> A4> 2 or 
/i 2 ~ A4> 2 , both Ap and r either one of them is not in agreement with the Planck and WMAP9 


experimental results. For example, taking A = 30mp, and random values of other parameters 
10 -10 mp < \y?\ < 10 4 mp, 10 -15 < |a| < 10 3 , 10~ 10 < |A| < 1 we get |A^J > 10 3 . If we assume 


that /i 2 -C A4> 2 , the potential (64) can be rewritten in simple form 


4 > 


F($) = A'($ 4 + a'$ 4 ln-) 


where 


a'- a 

A “4’ 


/ 

a = 


+ 72A 2 


(69) 


(70) 


167t 2 A 

The coupling constant A is determined to satisfy the constraint on A^, while as the predictions 
for n s ,r,a are given for fixed values of a', A. Fig. [I] shows the predicted values of n s , r and 
a for A = O.lmp (green), A = 30mp (red), A = 50mp (pink), and A = 500mp (blue) in the 
range of —10 3 < a' < 10 3 with the number of e-folds N = 60. We can see that for A = O.lmp and 
A = 500mp, r runs out of experimental region for almost values of a’ in the range —10 3 < a' < 10 3 . 
For A = 30mp, we need to require a' < —36 or a' > 6 to make sure n s and r are in agreement 
with experimental results (T6j, n s £ (0.94,0.98), and r £ (0.001,0.15). 

If we vary a!, A in the parameter region satisfying experimental results, the order of < > and 

the inflaton mass mostly does not change. From now on we take a' = — 10 2 ,A = 30mp for the 
below numerical calculations. 

From the minimal potential condition, we get < >~ 23.6 mp. The inflaton mass is calculated 
by the second derivative of the effective potential at the minimum. For Z N and VkM, the mass 
arises from (59) with notice that T = \f21Z[(p\. We obtain 


m$ 


= \fV"(<&) !$=<$> ~ 2.67 x 10 13 GeV, m Z N = 2 < $ >, 


m v iR = -y/^K < $ > . 


( 71 ) 
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FIG. 1: r vs. n s (upper panel) and a vs. n s (lower panel) for A = 0.1 mp (green dot), A = 30 mp (red 
minus), A = 50mp (pink plus), and A = 500mp (blue asterisk) in the range of —10 3 < a' < 10 3 , with the 
number of e-folds N = 60. 


Now let us calculate the reheating temperature. In this model, the inflaton couples to pair of 
Higgs, pair of gauge boson Z N and pair of Majorana neutrinos. We assume that < m Z N , 
hence the inflaton cannot decay into pair of Z N . The inflaton can decay into pair of Higgs with 
the decay rate 

A?!,.,,.,, < $ > 2 

r($ -> hh) = °’ ’ - . (72) 

v ’ 327rm$ v ’ 
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If the mass condition is allowed, the inflaton can decay into pair of vm 

r(® -> w iR ) = (73) 

If |Aio ; ii;i 2 1 3> —we get r(<f> -» hh) S> T(<h — > ViRym). The inflaton dominantly decays 
into pair of Higgs, therefore the reheating temperature is estimated as 

Tr = ^2 g* ^ (r$mp )2 ~ 10 20 GeV x |Aio ; ii ; i 2 1 , (74) 

where g* = 106.75 is the number of degrees of the freedom active at the temperature of the 
asymmetry production. The constraint T(3> —>• hh) <C requires |Aio ; ii ; i 2 1 <C ^ 10~ 6 . 

We find the limit Tp(max) < 10 14 GeV. Taking |Aio ; ii ; i 2 1 ~ 1CT 11 then Tr ~ 10 9 GeV satisfying 
the upper bound on reheating temperature to prevent gravitinos problem. In this case the thermal 
leptogenesis scenario may work to explain the baryon asymmetry. 

In other case, we assume 


therefore, 


|Cl < ——-~ 3.33 x 1(T 7 < |C 2 | ~ |dI, 

V2 < d* > " 


(75) 


m UlR <C m<s> < m U2R ~ rn„ 3R . (76) 

If |Ai 0; ii; 12 | < we get T(<P ->• hh) < T(<f> ->■ virvir). When Ai 0; ii;i 2 are negligibly 

small, the inflaton dominantly decays into pair of vir. The produced reheating temperature is 
given as 

Tr = ( f y) 4 ( r * m ^ ~ 1014 X Id- (77) 

This temperature is much lower than the RH neutrino mass since < > is at Planck value. 

We can apply non-thermal leptogenesis scenario, in which the v \r is produced through the direct 
non-thermal decay of the inflaton T. 


IV. LEPTOGENESIS IN THE 3-3-1-1 MODEL 

First, we consider the scalar sector. The scalar mass spectrum is considered in EE], in which 
—/~o;~As>u~u. In this work we assume the considered model contains 

• There are 9 Goldstone bosons A 4 , Gz, G' z , Gx, G* x , Gy, their interactions can be gauged 
away by a unitary transformation. 
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• One higgs gains mass at the electroweak breaking scale. This is the lightest massive Higgs 
bosons H and is identified as SM Higgs. 

• There are 9 new Higgs bosons namely, A, H\, H 2 , Hf, , H', H'*, which are heavy at the 

u> scale, while the mass of H 3 is proportional to A. 

In the gauge sector, let us collect the new gauge bosons beyond the SM. In the limit A>w> 
u, v, we get 

• One super heavy gauge boson ~ with the mass m 2 zN — Ag 2 N A 2 . 

• All the other new gauge bosons, Z 2 , X°, X°*, Y±, have mass in order 0(u). 

The lepton number of particles are considered in 019 ]. In particularly, the SM particles have a lep¬ 
ton number as usual. The new particles (Gx, H'*, , Gy . X° , Y~ ) have the lepton number equal 

to one, their complex conjugate have the lepton number equal to minus one while the remaining 
Higgs and gauge bosons have zero lepton number. 

Now in order to account for leptogenesis, we have to verify the lepton number violating inter¬ 
actions. Seeing that the lepton number L and baryon number B are conserved by VEVs of rj, y, p 
as mentioned in |8j. All interaction terms appeared after symmetry breaking in the considered 
model are conserved the lepton number. Hence it is clear that B, L violating number interactions 
should be broken in other way in order to explain neutrino mass and mixing as well as the matter- 
antimatter asymmetry of the Universe. The lepton number only can be violated in the interactions 
of Majorana neutrinos with non-zero lepton number particles. 

Let us remind the seesaw mechanism that explains the tiny neutrino mass and large mixing. 
The Lagrangian relevant to the neutrino mass has a form as 

Liz-mass — T T H.c). (78) 


The left handed neutrinos couple to the right handed neutrinos through the first term of the Eq. 


(78) and have a Dirac mass as 




(79) 


while the right handed neutrinos couples to themselves through the second term given in the Eq. 


(78) and have a Majorana mass as 


\m 


= -V2hUA. 


(80) 
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Hence, we can explain the smallness of the light neutrino masses via a type I seesaw mechanism 
|8] and predict six Majorana neutrinos as mass eigenstates, three heavy neutrinos iy t M and three 
light neutrinos I'iE, 

vm = ViR + i'I r\ rn UM = m™ = -y/2h'A, (81) 

ViE = v iL + vf L ; mf E = -m° (m") _1 (m^) T = 7 ^^h u {ti u )~ 1 (h u ) T . (82) 

We note that the considered model also contains three new neutral fermions N aR . They obtain 
the Majorana masses [8] via an effective interaction as 

J^{^aL)m{i’bL)n(XmXn)* T H.C. (83) 

The Majorana masses of the neutral fermions N aR are given 

[m NR } ab = (84) 

and the Majorana fermion states are 


Ni = N iR + N? R . 


(85) 


Based on the Majorana fermion states given in Eqs. 
Lagrangian C VR including the Yukawa terms in Eq. 
ViRi^D R v iR as follows 


(81), (82) and (85), we can 


(78) and the gauge-fermion 


rewrite the 
interaction 


£ 


V R 


{hab^aL’WbM + K b U aM P R U bM (j) ~ ^\jn^]ab^aMP r^M + H.c) 
+9N^iMl^ PRViM + ViM'iy l d ll P R V l M + H.C. 


( 86 ) 


To rely on Higgs physical states mentioned above, we obtain the physical interaction terms that 


violate the lepton number. In particularly the lepton violating interactions appeared in Eq. (86) 
are: eP R u M H b , NP R H'v M - 

We would like to emphasize that the lepton number violating terms also appear via the inter¬ 
actions of the light Majorana neutrinos, namely, eYglT - , NveX®* . However these interactions do 
not generate baryon asymmetry by m ■ In brief, this model contains the lepton number violating 
interactions, which are e P R vm , NP R H ' gy, cueW~ , NveX 0 * ■ We consider leptogenesis sce¬ 
nario at the temperature Tp satisfying Tp 2> ITeV. It implies that only vm can generate lepton 
asymmetry. 
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vertex 

coupling 

vertex 

coupling 

VaEVbMH 

p„ 

V2Au 2 +v 2 H 

VaEVbMH i 

p„ 

\/2\/« 2 +« 2 H 

e a v bM H^ 

"Kb _P„ 

\Zw 2 +i; 2 ^ 

NaVbAlH' 

U Kb p 
y/u 2 -\-UJ 2 ^ 

VaMVbMHs, 

h ^Pn + ^Pl 

V a M V b M Zj/ 

gNi^PR 


ufi lb p 

N a e b H+ 

p 1 VLj\ ab p 

\/u 2 +V 2 ” 

Vv 2 +u 2 R 1 V2jv 2 +u 2 M L 

N a v bE H ' 

uui\ ab P 

y/2s/u 2 +u 2 M L 

e a v b EW~ 

-X* 

e a e b H 

vh ab P„ _1_ D 

&a,e b Hi 

U ^ab p | p 

V2v'« 2 +t> 2 n/SvV’+u 2 ^ 

V2Vu 2 +v 2 R ‘ r N/SVii 2 -!-!; 2 L 


gswY 

e a e b Z£ 

~ gkAl 5 ),k = 1,2,N 

VaENbXn 


e a N b Y~ 


N a N b Zl 


NaN b Zj/ 

IgNi^PL 

N a N b H 2 

UJXab JD _i_ 0J X ba -p 

2 M L 2 M R 




TABLE I: Non-zero couplings of fermions appearing in loop diagram of vm —► e + and vm —► N + H'* . 


Before calculating the CP asymmetry of u 0 m , for convenience, we list all non-zero couplings of 
fermions appearing in loop diagram of vm ->e + and um —> N + H'*. 

All possible one-loop diagrams, which can contribute to the CP asymmetry from the decay 
vm —> e + H§ are listed in Fig. [ 2 J The interference of the tree level and one-loop level (2a, 2b), 
(3b) with the propagator ujm, (6) with the propagator (H^,ei) gives dominated contribution to 
the CP asymmetry. We obtain 


j(l) _ r {vkM ~> + #5~) — r {vkM ~> Cj + H 5 ) 


" v kM 


2r 


v kM 


87rCo 


r ) 9^w+H 1 H~ + eX AH+H7 + ( 9lV + 91a)^ Zi H+H~ + ( 92V + 92a)^z 2 H+H~ 


+ ( 9NV + 9na)X Z nH + h - (4 - 2g z log[\ + l/g z ]) s} E Im [h^ 


E t 1 - (! + 9j)log[l + 1/9:} + (1 - 9:)~ l ] (87) 


where T UkM is the total decay rate of VkM at tree level, 


9z = 9j = 


m 


2 ’ 
v kM 


m 


'W+H.H- - 2 > "AH*Hi 


2 ’ 
VkM 


Co — (2 + s 2 p) ^2 I h/k 1 2 — (2 + Sp)(h u lh u )kk, tp — v/u, 


'ZiH+HI 


gc2w 

2 c-w 


X z 2 h+h- - 


g{c 2 0 - C2WSp) 
2 cw ^3 — 


A 


ZnH+H- 


9N 2/3 , ( ffiy + 5lj4 ) = -(g 2 y +52A) = (dNV + 9na) = pp- (88) 

3 2cp/ 3 


3 


3 
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Here we ignore the mixing between Z * and Z 2 since w >> u, v. 

Now we consider CP asymmetry of the decay I'm —> N + H *. All possible loop diagrams are 
listed in the Fig. [3j The interference of the tree level and one-loop level (2c), (3) with the propagator 
UjM, (6) with the propagator (H b ,ei) gives dominated contribution to the CP asymmetry. We 
obtain 

i(2) _ r (vkM — t Ni + H'*) — V{y kM —> Ni + H') 




2 T 


v kM 


1 


87tCo 


g = A z 2 h'*h’ + ^9n^z n h'*h'(^ - 2g z log[l + 1 /g z }) 

,2 3 


3 — 4s 


w 


^ 2 lm[h^h 


i 


+ g 7r( o X/ [l _ (1 + 9j)l°9[ 1 + 1/ffj] + s |(l — 9j) *] I m[(h^h u ) k jh" k hy\, (89) 

j 


where 




H'*H' = 


gcw 


3 4s^ 


Az„ 


H'*H' = 




(90) 


We would like to notice that since the coupling A e aUbM H~ = s ^ab^ > R while ^ aUbM jjt — h v ab PR , 
the factors s|, s^ appear in (87) while s°, s 2 appear in (89). In this work we take u ~ u and thus 
S/3 = 1/V2. 


Let us comment on neutrino mass and mixing. The light neutrino mass matrix is given Eq. (82). 
In order to diagonal this matrix, we have to use the U matrix. It is nice to note that the lepton 
mixing matrix was studied by the Pontecorvo-Maki-Nakagawa-Sakata (PMNS). The standard form 
of this mixing matrix is given 


/ 


CpMNS = 


C12C13 


S12C13 


si3e 


— i5 ^ 


iS 


-312^23 “ Ci2S23Sl3e C 12 C 23 — Si2S23Sise $23^13 


(91) 


\ S12S23 — Cl2C23'Sl3e*' S ~C12S23 — Sl2C23'Sl3e l<5 C23C13 


where Cij = cos 0 jj , Sij = sin Qij and the values of 0^ are determined by the global analysis [18] , 
namely 

sin 2 023 = 0.466^;SSS; sin2 = 0 - 312 -aoo?8 0 o°049; sin 2 = 0-016 ± 0.010(< 0.046). 

(92) 

5 is unknown CP violating Dirac phase. 

On the other hand, the square of charged lepton mass matrix and light neutrino mass matrix 
are diagonalized by two unitary transformations 


MiUi = Diag(m 2 , m 2 m 2 ); t/J m e J^U v = Diag(m vi , m V2 ^ m 


’’VZ/ 


(93) 
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The C/pMNS is defined as 


^PMNS-P = U} U u , 


(94) 


where 


P = 


1 0 0 
0 e ia 0 
0 0 e ip 


(95) 


where p , a are CP violating Majorana phases. If we ignore the mixing between the charged lepton, 
then we can get 


U v = CpMNS-P- 


(96) 


We assume that h' v = Diag^h'^, h'^, h. 33 ) then m^ 1 = Diag(m l , 1M ,m U2M ,m l , 3M ). Using the analysis 
in [IS], the most general h u matrix is given by 


h u = — Diag(^/m UlM ,y/m U2M , y/m„ 3M ).R.Diag(y/m^, s/rn^, 


(97) 


where R is orthogonal matrix expressed in terms of arbitrary complex angles 6\ , 62 , 63 as following 


( C2C3 -ClS3 - S1S2C3 SlS3 - C1S2C3 \ 


R = 


C2S3 C1C3 — S1S2S3 —S1C3 — C1S2S3 
\ S2 SlC‘2 CiC 2 

where c % = cos 6 i, s* = sin d,, i = 1 , 2 ,3. 


(98) 


From the Eq. (97) h v 'hy has the form 


h^h u = ~^U.Diag( y /rfhP 1 , y/m v2 , y/m U3 ).. Diag(m UlM ,m U2M , m„ 3M ).R.Diag{y/m Ul , % /m^').t/ t . 

(99) 

For the light neutrinos masses, we fit the experimental results 


Am^ = rnP - ml = 7.53 x 10 5 eV 2 , A rrtL = = 2.44 x 10 d eVC 


-3„w2 


( 100 ) 


The asymmetry now can be considered as function of the phase S,p,a , the heavy 

majorana neutrinos masses and the complex angles 01 , 02 , 63 . For simplicity, we assume 
6 \ = 62 = #3 = 0- In this work we consider the CP asymmetry due to the decays of the lightest 
heavy Majorana v\m- The detail will be presented in the subsections below. 
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H,Hi ' i H+ 




Hiy | ni Wjt ' I 










(x, y = (H/H i, v lE ), (//', JVj), (H 3 ,vim), {Z*,v lM ), (H“ , e,)) 


FIG. 2: One-loop diagram contributing to the asymmetry from the decay vm —t e + 


The baryon asymmetry and lepton asymmetry are given as 

n B ~ ng 


VB = 
V L = 


s 

Hi - Tlj 


( 101 ) 


where s is the entropy density. The lepton asymmetry can be transformed into a baryon asymmetry 
by non-perturbative B + L violating (sphaleron) processes [20], giving 


VB = a(vB ~ Vl) = 


a — 1 


VL, 


where 


a = 


8 n g + 4 riH 
22 n„ + 13 tih ’ 


( 102 ) 


(103) 


with riH is the number of Higgs and n g is the number of fermion generations. We get 


VB = 


8 n g + Ann 
14 n„ + 9 nu 


VL- 


(104) 
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H' 



H 3 ' i H' 


1h> 



N< 


(Id) 



Hi 


-T~ 

I 


H' 





+^4 + 


Ni 


(lb) 




H' 






S N 


VkM V y VjE 

ei(N t ) 

(7) 

(x ,y = (H/H^vie), (//', IV,), (H 3 ,v lM ), (Z™ ,v lM ), (H“,e,)) 


* 


Ni 


FIG. 3: One-loop diagram contributing to the asymmetry from the decay vm —> N + H'*. 


As the analysis in |2T] , taking n# = 2 and n g = 3, we get 

8 , , 
11B = L- (105) 

Now let us calculate tjl in thermal and non-thermal leptogenesis scenario. 

A. Thermal production 

In the thermal scenario, the heavy Majorana neutrinos are produced in a thermal bath. At 
T > m UlM , the CP asymmetry generated by v\m decays can be washed out due to inverse decays 
and scattering processes. That why the CP asymmetry is weighted by the washout efficiency. 

For the channel ■> ^ ie CP asymmetry depends on flavor because Li{ei) = 1. 

However, since L(Ay) = 0, L(H') = —1, the CP asymmetry due to the decay UkM NiH'*, NiH' is 
considered flavor independent. The Boltzmann equations for the lepton asymmetry can be divided 
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2(1) 

by two forms, one is the equation for the flavored lepton asymmetry corresponding to and 

another is treated by the conventional computation for e$ m = • 

The interference of the tree level with loop diagrams contained gauge propagator is vanished in 
summation of all the indexes i,l = 1,2,3 if the CP asymmetry has the same weight for all flavors, 


E M*K] = ME'«] = o- 


(106) 


Therefore, from Eq. (89) 


e ( 2 ) = 
VlM 


E 


r®(2) . 


1 


'VlM 


8itCo 


V9j [! - (! + 9j)log[ 1 + 1/gj] + s 2 ( 1 - gj) *] lm[{(h^h u ) lj )‘ 2 


1.6 x 10 -2 ^ 


y/gjQrfhy) n 


(107) 


Eq. (87) can be reduced by taking g = 0.65, sfy = 0.231 as 


=-*(!) = 
C t / l M 


-1.6 x lO" 4 E; Im [hfihtt - 0.6 x lO" 2 V- Jg 






(h^h v ) 11 

In the thermal leptogenesis, the washout parameters are defined as 

VjyiM -* ejH+, ejHj[) = s}h^h^m UlM 
H(T = m VlM ) 8vr 

T(^im) _ ( 2 + s 2 p)(h^h l ')u m V\M 


(108) 


Ki = 
I< = 


H(T = m UlM 

By varying 5, cr, p, 8 for m UlM ~ 10 9 GeV, m 


8 tt 


(109) 


"2M 


m 


V3M 


10 3 m UlM we figure out K 3> 1 


for all values of CP parameters 5, a, p, and the complex angle 8 . The lepton asymmetry can be 
approximated as given in [22]. In the strong washout regime, K 1 


rp L ~ 0.3 


-( 2 ) 

'UM 


g* 


0.55 


Vl 


J(i) 

-um 

g* 


K 
8.25 


1.16 


+ 


I An | Ki 

0.2 


1.16' 


-1 


( 110 ) 


where A n = -151/179, A 22 = A 33 = -344/537. 

The baryon asymmetry is related to the lepton asymmetry as 

8 


VB = 


15 


( 8l + Vl)- 

i= 1,2,3 


(111) 


From all expressions above, we see that ijb depends on 5,cr,p, and 8 . The baryon asymmetry gs 
in the region (5 x 10 —11 ,10 _1 °) on the plan of the complex angel 8 is shown in Fig. [ 4 ] for 


5 = 4.3 rad, a = —1.5 rad, p = — 1 rad, 

m v2M = m u3M = 10 3 m UlM , m UlM = 10 9 GeV, m Ul = 0.01 eV. (112) 
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The red regions indicate that in order to satisfy 5 x 10” 11 < t/b < ICR 10 , we need to require 
— 1.01 < Im[0] < 1.8 when varying Re[0]. The limit of lm[0] keeps the same if we extend the range 
of Re[0]. The rjs is considered as function of pure imaginary 6 (red) and pure real 0 (blue) as 
shown in Fig. [Hj We see that r]B changes a lot when varying pure Im[0] while it seems to keep the 
same order when the pure Re[0] alters. The baryon asymmetry varies little as function of the CP 
phases presented in Fig. [6j 

If we study the case m UlM = 10 9 GeV,m U2M = = 10 5 m UlM , the constraint on the complex 

angle 6 is stricter in order to satisfy the experimental results on baryon asymmetry. 



R e[0] 

FIG. 4: Contour plot of rjs in the region 5 x 10~ n < ijb < 10~ 10 on the plan of the complex angel 9 for 
<5 = 4.3 rad, cr = —1.5 rad, p = — 1 rad, m V2M = m„ 3 M = 10 3 m UlM , m VlM = 10 9 GeV, m Vl = 0.01 eV. 


B. Non-thermal production 


In the non-thermal scenario the reheating temperature can be lower than the lightest heavy 
Majorana. The total CP asymmetry is the summation of all flavor CP asymmetry, 




B 


r»(l) 


£ VkM e VkM 


m \ - 


(h^h v 


Ikk 


(113) 
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e e 


FIG. 5: tjb vs. pure imaginary 9 (red) and pure real 9 (blue) for S = 4.3 rad, a = —1.5 rad, p = — 1 rad, 
m„ 2M = m„ 3M = 10 3 m VlM , m VlM = 10 9 GeV, m Vl = 0.01 eV. 



S (rad) 


cr (rad) 


FIG. 6: pb vs. S (left) and tjb vs. a = p (right) for 9 = 0.87 1 (red) and 9 = —0.18/ (blue), and other 


parameters given in (112). 


where 


Bj = 


1 


87 t (2 + s 2 p) 

11 

16 0 TTy/g] 


Vs] [( s /3 + 1)(1 _ (1 + + 1 / 9j]) + s /3( s /3 + 1)(1 — 9j) 


-n 


(114) 


The lepton asymmetry is related with the CP asymmetry through 


3 T r 

Vl = 7 ze VkM x Br k x ——, 

2 ?TlcJ> 


(115) 


where Br^ denotes the branching ratio of the decay channel <f> —> I'kM^kM- 
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As analysis in the previous section, we assumed that rn V] M <C < m U2M ~ m^ 3M , < m Z N 

and T(<h —> hh ) -C r(4> — > i'irRir) when Aio ; ii ; i 2 are negligibly small, therefore, 


VL - -£ 


vim 


X 


Ik 


Combining Eqs. (77, 105, 113) with 'Uj U r^j 174 GeV we get 

m- 
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b V\ M 


T]B — 0.4 X 
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j= 2,3 m v 


\/2 < <h > 
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'jM 
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(117) 


with notice that the formula h^h u given in Eq. (99). Putting 


5 = 4.3 rad, a = —1.5 rad, p = — 1 rad, 6 = 1.46/, 

m„ 2M = 3M = 10 3 m UlM , m UlM = 2.34 x 10 11 GeV, m Ul = 0.01 eV, 

= 2.67 x 10 13 GeV, < $ >= 23.6m P , (118) 


we get 


Pb — 8-92 x 10 


-ii 


(119) 


This value of baryon asymmetry is in agreement with [|23j, tjr = (8.75 ± 0.23) x 10 -11 . 

Let us consider how ijr depends on the complex angles and CP phases one by one. Fig. [7] 
shows t]r in the region (5 x 10 —11 ,10~ 10 ) on the plan of the complex angel 6 for m^ 2M = m,„ 3M = 
lO 3 ??!^^, m UlM = 10 11 GeV (red), and m„ 2M = m„ 3M = 10 5 mi, 1M , m UlM = 10 9 GeV (blue) and 


all other parameters as given in (118). We see that in the red region —2.05 < Im[$] < —1.68 or 
1.49 < Im[0] < 2.28 and in the blue region lm[0] ~ 3.3 or lm[0] ~ —3.4 when varying Re[0] even 
though if we extend the plot range for both axes. It means that it is free to choose the value of 
Re[0] but lm[0] is quite a strict constraint. rjB depends strongly on lnr[0], while it changes lightly 
when varying Re[0]. This conclusion is more clearly in Fig. [8j in which rjB is considered as a 
function of pure imaginary (red) and pure real (blue) 8 . 

Fig. [9] shows t]b as a function of Dirac CP phase 5 (left) and Majorana CP phase cr = p (right) 


for 9 = 1.46/ (red) and 9 = —1.46/ (blue), and the choice of other parameters given in (118). In 


brief, we see that ijb does not depend much on the CP phase but depend on the imaginary of the 
complex angle 9. This conclusion is the same as analysis in thermal scenario. 


V. CONCLUSIONS 


We have studied generation of inflation and leptogenesis in the 3-3-1-1 model by considering the 
symmetry breaking of the XJ(1)n gauge group at the GUT scale. The model contains two super 
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FIG. 7: Contour plot of ijb in the region (5 x 10 11 < rjs < 10 10 ) on the plan of the complex angel 6 for 
rrii ,im = 10 11 GeV (red) and m VlM = 10 9 GeV (blue). 



FIG. 8 : r/B vs. pure imaginary (red) and pure real (blue) 9. 


heavy particles with mass proportional to A, the new gauge boson Z N embedded to U(1)n and the 
scalar Higgs boson H% ~ S 4 . All other new massive particles get mass in order of u>. The singlet 
Higgs <j) with < > at the GUT scale can play the role of inflaton. The quantum corrections to 

the potential of inflaton is taken into account, thus there appears logarithm function of inflaton, 
making the presently considered model’s inflation different from chaotic one. In this work, we have 
figured out the parameter spaces appeared in the inflaton potential matching the experiment on 
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S (rad) 


cr (rad) 


FIG. 9: r) B vs. S (left) and t/b vs. a = p (right) for 0 = 1.46/ (red) and 9 = —1.46/ (blue). 


the spectrum index n s , the tensor to scalar ratio r, the running index a as well as the amplitude 
of the curvature perturbation A^. The inflaton mass is obtained in an order of 10 13 GeV. 

After the inflation, the heavy Majorana can be produced in a thermal bath or by decay of the 
inflaton. Depending on the Higgs couplings Aio,11,12 in comparison with the Yukawa couplings h 
leptogenesis is considered in thermal or non thermal scenario. We have shown how the 3-3-1-1 model 
generates lepton asymmetry then converts into baryon asymmetry in both cases. It is interesting 
that the model contains an extra channel contributing to the CP asymmetry. The heavy Majorana 
particles can decay into neutral neutrinos A r , and neutral complex Higgs H' with the coupling 
different by factor sp from the original channel, i>kM —> efH^. In thermal leptogenesis, the CP 
asymmetry generated by the new channel is considered flavor independent, while the ordinary 
channel is treated as flavor dependent due to the different lepton number of N t and e*. It leads the 
interference of the tree level with loop diagrams appeared gauge propagator to contribute to the 
CP asymmetry for the decay —> efH^. This feature is new compared to other leptogenesis 

models. 

The thermal and non thermal leptogenesis have been calculated in detail. In order to get non 
zero CP asymmetry we need to consider the complex Yukwa coupling matrix h u by expressing it 
in terms of the neutrino mass and mixing matrix and the orthogonal matrix R. We have presented 
how tjb depends on the CP phases 5 } a,p and complex angle 6 = 6 \ = 62 = #3. The baryon 
asymmetry is not much sensitive to the value of CP phases or pure real 6 but it alters a lot as a 
function of pure imaginary 6 . This property is the same for both leptogenesis scenarios. Thank to 
the orthogonal matrix R and the complex angle 9, which makes the baryon symmetry completely 
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in agreement with the experiment for both cases. One different thing of the two scenarios is that 
at any point of Re[0] we always can find lm[0] satisfying a fixed value of tjb in non thermal case, 
but there is restriction of choosing pair of (lm[0], Re[#]) in thermal scenario to match experiment 
on t]b. We know that the baryon asymmetry depends much on lm[0] and it is easy to see that 
from the Fig. [5j there is an upper limit on the baryon asymmetry if we consider tjb as a function 
of Im[$] in thermal scenario because of the effect of washout efficiency. However, there is no upper 
bound for 77 ^(Im[0]) in non thermal case, see Fig. [ 8 j By considering non thermal leptogenesis, 
the reheating temperature Tr can be reduced much lower than the lightest heavy Majorana mass. 
In brief, the 3-3-1-1 model at the GUT scale successfully explains the baryon asymmetry of the 
universe by studying both thermal and non thermal leptogenesis mechanisms. 
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